#---------------------------------#
# Bayes model fits          #
#---------------------------------#
rm(list = ls())
load("model_output/multi_factor_model.RData")

plot(multi_factor_model)